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^ ; Abstract 
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A supersymmetric inversion method is applied to the singlet 1 Sq and 1 Pi 
neutron-proton elastic phase shifts. The resulting central potential has a one- 
, pion-exchange (OPE) long-range behavior and a parity-independent short- 

range part; it fits inverted data well. Adding a regularized OPE tensor term 
also allows the reproduction of the triplet 3 Po, 3 Pi and 3 Si phase shifts as well 
^ ■ as of the deuteron binding energy. The potential is thus also spin-independent 

(except for the OPE part) and contains no spin-orbit term. These important 
, simplifications of the neutron-proton interaction are shown to be possible 

O ! only if the potential possesses Pauli forbidden bound states, as proposed in 

' the Moscow nucleon-nucleon model. 
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Although of fundamental importance for the whole field of nuclear physics and despite 70 
years of intensive research, a definitive form for the nucleon-nucleon interaction has not yet 
been established. Since an individual nucleon is made up of three quarks, this interaction 
should in principle be deduced from an ah initio quantum chromodynamics description of 
a six-quark system. However, such a six-body highly-relativistic problem is out of reach of 
our present calculation capacities. Hence, phenomenological interactions, i.e., interactions 
constructed to reproduce nucleon-nucleon experimental data (deuteron bound state and 
scattering cross sections) have to be constructed. These interactions may be divided in two 
categories. 

The first category is based on the theoretical framework of meson exchange and fulfills 
the symmetries of the system. The best known are the Paris jl], Nijmegen 0, Bonn [Q, and 
Argonne j|] interactions. All these models agree on the long-range part of the interaction 
(r > 3 fm, where r is the inter-nucleon distance): it is well described by a central (C) and a 
tensor (T) one-pion-exchange (OPE) potential. Both terms depend on the total spin S and 
parity tc of the nucleon-nucleon system. The whole potential reads 

V = V s - r - + Vg PE (S, tt) + V^ PE (S, tt), (1) 

where the short-range part ^ s r ' follows from the exchange of other mesons and from a purely 
phenomenological contribution. Its precise form differs from model to model but it at least 
contains a central, a tensor, a spin-orbit (SO), and a quadratic spin-orbit (QSO) term, 

V**- = V^iE c . m .,S,ir) + V*- T iS,ir) 

+ViSiS,n) + V^ (S,K), (2) 

all depending on spin and parity. The central-part dependence on the center-of-mass energy 
-E c . m . makes the potential non local. 

The second category of interactions is physically less justified: the long-range part is 
also described by the OPE potential while the short-range part has a simpler but more 
phenomenological structure. In the Reid potential 0, V s ' v ' has neither QSO term nor 
energy dependence. The price to pay for this simplification is a partial-wave dependence of 
the potential: instead of a dependence on S and tt, it depends on S and on the orbital and 
total angular momenta L and J. There is thus an implicit non-locality and this interaction 
does not fulfill the symmetries of the system. An even simpler form is proposed in the latest 
version of the Moscow model ||. In this case, the short-range potential consists of a local 
part only depending on spin and parity, without tensor term, 

V^ l = V^iS,rr) + V^iS,rc), (3) 

and a non-local separable term depending on the partial wave V££(S, L, J). This further 
simplification of the interaction structure is made possible by the introduction of Pauli 
forbidden bound states (PFSs) in the lowest partial waves. These non-physical states are 
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aiming to simulate the Pauli exclusion principle between the six quarks constituting the 
interacting nucleons [0]. Their presence results in a deep potential without repulsive core, 
whereas other nucleon-nucleon interactions, which only possess one physical bound state 
(the deuteron), are shallow and present a repulsive core below about 1 fm. 

In the present work, we show that the use of a deep potential allows a further simpli- 
fication of the short-range local part (H). Neither the spin and parity dependence nor the 
spin-orbit term seem to be necessary to fit data, so that one only has 

V£ c r al = V^. (4) 

Spin and parity dependences, or equivalently spin-spin and isospin-isospin terms, are gener- 
ally considered as essential ingredients of the nucleon-nucleon interaction, in contrast with 
Eq. (^). Besides the presence of PFSs, the present simplification requires a more complicated 
central part than in the Moscow model, where an exponential or Gaussian form factor is as- 
sumed from the beginning. Here we deduce this form factor from neutron-proton scattering 
data through the use of the supersymmetric inversion method of Ref. ||. This technique 
allows the construction of a partial-wave-independent potential from elastic scattering phase 
shifts of several partial waves. In Ref. || , the inversion method is applied to the system of 
two a particles, for which the inclusion of PFSs has been known to considerably simplify 
the interaction for a long time. We show here that a similar simplification occurs in the 
interaction between two nucleons, a fact ignored up to now. 

Our inversion method is a two-step procedure ||. We first invert the phase shifts 5(k) of 
the sine; let (S = 0) 1 Pi partial wave. As "experimental" data we use the multienergy phase 
shifts PWA93 of Ref. g] at E c _ m , = 0.05, 0.5, 2.5, 5, 10, 15, ... , 140 MeV (below the inelastic 
pion-production threshold). The wave number k (in fm -1 ) is related to -E c . m . (in MeV) 

through k = \j2fiE c _ m Jh 2 , where fi = 469.45926 MeV/c 2 is the neutron-proton reduced 
mass with he = 197.327053 MeV fm. We use a rational approximation of the scattering 
S'-matrix 

S(k) = e 2 ^ « S° PE (k) J} (5) 

771=1 Km ™ 

where Sq PE is the S- matrix of Vq PE (V® PE vanishes for singlet waves) and K m are complex 
wave numbers fitted to experimental data. Potential Vq PE depends on spin and parity 
through the spin and isospin Pauli matrices of the nucleons cr 1 , <x 2 , Ti, t 2 . It reads (in 
MeV) 

V° PE (S,rc,r) = 0.075 m^c 2 $ c (r) (<r t ■ <x 2 ), (6) 

where an averaged pion mass m n = 138.0363 MeV/c 2 has been used (distinguishing between 
neutral- and charged-pion exchange does not qualitatively modify the present results). A 
dipole-regularized form factor $ c is used [Eq. (9) of Ref. 0] with a cutoff parameter Ac 
given below. Other form-factor choices (as in Ref. ||) will be discussed elsewhere. For the 
1 Pi partial wave, one has Ti • r 2 = (T\ ■ <r 2 = —3 in Eq. (H). The corresponding OPE phase 
shifts are represented by a dotted line in Fig. [1]. They are satisfactory at low energy only. 
With three additional S'-matrix poles [M = 3 in Eq. ([5])] located at K m = ±1.848 + 1.860? 
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and —1.848? fm 1 , a perfect data fit is obtained on the whole energy range (upper dashed 
line in Fig. |l|). The values of the poles have been calculated with the method presented in 
Ref. |g . 

The potential corresponding to the rational S'-matrix (|5|) is a generalized Bargmann 
potential ||11|| . It may be constructed by three successive super symmetric transformations 
|i~2|l of the l Pi radial Schrodinger equation with the OPE potential. It has no bound state 
and is purely repulsive (upper dashed line in Fig. ^) in agreement with traditional shallow 
nucleon-nucleon interactions. To calculate 1 5'o phase shifts, we modify the OPE term with 
T\ • r 2 = 1 instead of —3 in Eq. (g) and we keep the same short-range part [Eqs. ([]]) and 
(||)]. This makes the central potential slightly attractive above 1.5 fm (lower dashed line in 
Fig. 0). However, the repulsive Vq 1 ' dominates and the l So phase shifts are negative (lower 
dashed line in Fig. [I]), in disagreement with experimental data. Modifying Ac does not solve 
the problem. This explains why shallow nucleon-nucleon potentials always have a strong 
parity dependence: a shallow potential is unable to fit both the 1 S and 1 Pi phase shifts 
simultaneously. 

Let us now perform the second step of our inversion technique [§], namely, the addition 
of a Pauli forbidden bound state to the 1 P\ effective potential just obtained, in order to 
construct a deep potential. Such an addition may be done without modifying the 1 Pi 
phase shifts (phase-equivalent addition) with the help of two successive supersymmetric 
transformations [L| . Both steps of the inversion method are thus performed with the same 
algebraic formalism of supersymmetric transformations. Let us define the mean radius rpig 11 
of the added PFS by 

oo I 

\^ FS (r)\ 2 dr = -, (7) 

mean / 
PFS 

where ^pfs is the PFS normalized wave function. An important feature of the phase- 
equivalent PFS addition is that both the binding energy -Epfs an d rppg n may be chosen 
arbitrarily. Here we choose them to fit the 1 Sq phase shifts together with the 1 P\ ones. 
However, the obtained central potential has a singular r~ 2 attractive core at the origin and 
cannot be used as such to estimate 1 S phase shifts. Hence we regularize it as V reg + ar 2 + br 3 
below a regularization radius r reg , as in Ref. ||. This only slightly affects the X P\ phase shifts 
provided r reg is lower than the classical turning point for this partial wave at the highest 
considered energy. 

For E PFS = -101.9 MeV, rf%f = 0.6064 fm, V Teg = -2071 MeV, r reg = 0.8231 fm and 
A c = 468.1 MeV/hc, we get a very good fit of the 1 S and 1 Pi phase shifts (solid lines in 
Fig. [I]), with root-mean-square (rms) deviations of 3.9 and 2.3 x 10~ 3 7r respectively . The 
corresponding central potentials are deep and have no repulsive core. They are represented 
by solid lines in Fig. |2|, together with the energy and mean radius of their PFSs. Let us 
recall that, by construction, the s.r. term of the potential is parity independent since it fits 
both an odd ( x Pi) and an even ( 1 S'o) wave. Thanks to the soft cutoff parameter, the OPE 
term only introduces a small parity dependence. Such a soft cutoff is shown to be very 
satisfactory in Ref. || and seems to be a constant feature of deep potentials. 

In Fig. |3|, we compare the form factor of our potential (solid line) with that of the local 
part of the Moscow potential || (dashed line) for the 1 S'o partial wave. Both potentials have 
an OPE tail (dotted line) but have different short-range behaviors. In the Moscow model, an 
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exponential form factor is used, while the inversion potential has a more complicated form 
factor, which could not have been easily guessed from phase shifts. The inversion potential 
has no compact analytical expression but it is reasonably well approximated by a sum of 
two Gaussians 

V^-(r) = -3050 e- 5 ' 890r2 - 70.39 e" ' 5687 ^ (8) 

(in MeV, r in fm), as shown by the dash-dotted line in Fig. [| Adding an OPE term with 
Ac = 406.2 MeV /he to potential (|8|) provides rms deviations of 9.3 and 2.4 x 10~ 3 7r for the 
1 S'o and 1 P 1 phase shifts respectively, which is still in good agreement with experimental 
data. In the following, we only use the numerical inversion potential but we have checked 
that our conclusions also hold for the analytical potential which is of easier use for 
practical applications. 

Let us now show that the inversion potential is also able to fit triplet (S = 1) partial 
waves, in particular the deuteron bound state in the 3 Si- 3 Di coupled waves. We have found 
that it is sufficient to add a tensor OPE term to the above central potential to reproduce 
the experimental deuteron binding energy = —2.224575(9) MeV as well as the 3 Si phase 
shifts (Fig. [D with an rms deviation of 6.0 x 10 _3 7r. A dipole-regularized form factor $^ 
[Eq. (12) of Ref. 0] is used with a smooth cutoff parameter Af en = 601.14 MeV /he. The 
deuteron wave function is then similar to that of the Moscow potential (Fig. 5 of Ref. ||) and 
has a node in the S wave. The triplet-odd 3 Pq and 3 Pi phase shifts are also reproduced (Fig. 
f|, rms deviations of 1.3 and 0.9 x 10~ 3 7r respectively) with a cutoff parameter A^ dd = 517.7 
MeV /he. Finally, the mixing parameters have a good qualitative behavior but ei, €2 and €3 
are too small; this indicates that the tensor form factor could be slightly improved. 

The parity-independent central potential deduced from inversion of singlet-even and 
singlet-odd phase shifts is thus also valid for triplet-even and triplet-odd partial waves. 
Hence in addition to the parity, it is independent of the spin. To our knowledge, the present 
potential is the simplest phenomenological model of the nucleon-nucleon interaction able to 
fit all the above phase shifts simultaneously. Let us insist on the fact that the potential is 
purely local (no energy dependence, no J or L dependence) and has no spin-orbit term. The 
potential of Ref. || , on the contrary, has a spin-orbit and a non-local term for the 3 Pq and 
3 Pi waves. Moreover, it has no PFS for these waves, whereas our potential has one PFS for 
all S and P waves. The present model has thus a simpler structure. 

However, our potential does not reproduce the phase shifts of every higher partial waves 
(J > 2), which means that the present work is only a first step toward a simpler nucleon- 
nucleon interaction. In particular, we get negative 3 P2 phase shifts whereas the Moscow 
model agrees with the positive experimental values. Moreover, the D and F phase shifts 
have a satisfactory low-energy behavior but are too large at higher energies, as illustrated 
for instance in Fig. || (solid lines) for the singlet partial waves. A similar result is obtained 
for triplet waves. The higher partial waves (G, H, . . . ) are well reproduced because they 
mainly depend on the OPE tail of the potential. The discrepancy for the D and F waves is 
also encountered with the local part (|3J) of the Moscow potential (dashed lines in Fig. ^). In 
Ref. H, this is solved by the introduction of a repulsive pseudo-potential term V££ for these 
partial waves. More conventional non-local terms, such as an energy- dependent central part, 
could also be used and the possibility of a spin- and parity-independent non-local term could 
be considered. 
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As a conclusion, this work gives some hope for an important simplification of the phe- 
nomenological nucleon-nucleon interaction. A satisfactory local part has been found, which 
consists of a deep short-range central potential, independent of both spin and parity, and a 
regularized OPE potential (central and tensor terms). A spin-orbit term does not seem to be 
necessary, while the question of a satisfactory non-local term is still open. We plan to apply 
the same technique to proton-proton data (the Coulomb interaction does not complicate the 
inversion method |TTJ) in order to study the charge dependence of our model. We would 
also like to generalize the present inversion technique to the coupled-channel case, which 
would allow a deduction of the tensor form factor from experimental data. First steps in 



this direction have already been made |fH],|15 
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FIG. 1. Neutron-proton 1 Sq and 1 P\ phase shifts: data from Ref. Q (crosses), shallow- 
cashed lines) and deep-potential (solid lines) calculations. The l P\ phase shifts of both potentials 
are almost indistinguishable. The 1 P± phase shifts of the OPE potential (jGj) for Ac = 468.1 MeV /he 
(dotted line) are also shown. 
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FIG. 2. Shallow (dashed lines) and deep (solid lines) central potentials obtained by super- 
symmetric inversion for the 1 Sq and 1 P± partial waves. The 1 So and 1 P± potentials are almost 
indistinguishable in the left panel, where the energy (horizontal lines) and mean radius (vertical 
lines) of the deep-potential PFSs are also shown for both partial waves. 
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FIG. 3. Central potentials for the 1 Sq partial wave: dipole-regularized OPE potential @ with 
Ac = 468.1 MeV/hc (dotted line), deep inversion potential (solid line), deep Moscow potential of 
Ref. Q (dashed line), sum of the analytical s.r. potential (g) and of the OPE potential @ with 
A c = 406.2 MeV /he (dash-dotted line). 
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FIG. 4. Neutron-proton 3 Si, 3 Po and 3 Pi phase shifts: data from Ref. || (crosses) and 
calculations (solid lines) with the central inversion potential and an OPE tensor term. 
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FIG. 5. Neutron-proton singlet phase shifts: data from Ref. § (crosses), calculations with 
the inversion potential (solid lines) and with the local part @ of the potential of Ref. Q (dashed 
lines). 
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